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Abstract. The status of lattice calculations in the light hadron sector is 
reviewed. Special emphasis is given to recent lattice determinations of the mass of 
the strange quark. The impact of non-perturbative renormalization and control 
over lattice artefacts on the attainable precision is discussed in detail. Furthermore 
the influence of dynamical quark effects is assessed. 



1. Introduction 

Quantum Chromodynamics (QCD) is widely accepted as the theory describing 
the strong interaction. It is formulated in terms of quarks and gluons and 
contains relatively few free parameters, namely the gauge coupling g and the masses 
mu, m^, nis, ... of the various quarks. These parameters must be fixed by using 
experimental input: a typical example is the determination of the sum of quark masses 
(ttt-u + nig), using the experimental value of the kaon mass, mK+: 

mK+ ^ (mu + rus). (1) 

Clearly, in order to connect the two sides of this relation one has to "solve QCD" for 
the case at hand. 

The QCD couphng constant g, which describes the coupling strength, depends 
on the momentum transfer q between quarks and gluons. The well-known property 
of "asymptotic freedom" implies that g decreases for large q. By contrast, the 
coupling becomes large when momentum transfers of the order of typical hadronic 
scales are considered, say, q « 1 GeV/c. As a consequence, perturbation theory in g"^ 
becomes an inadequate tool to study QCD at low energies. Hence, a non-perturbative 
treatment is required in order to understand hadronic properties on a quantitative 
level. The formulation of QCD on a (Euclidean) space-time lattice provides such a 
non-perturbative method to compute the relations between experimentally accessible 
quantities and the parameters of QCD through numerical simulations. 

Here I restrict myself to QCD at zero temperature, focusing on the spectrum of 
light hadrons and the determination of the strange quark mass from first principles. 
A more general overview of the status of lattice calculations can be obtained from the 
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proceedings of the annual Conference on Lattice Field Theory and other review 
talks |-|. 

In the next section I shall outline the basic concepts of the lattice formulation. 
Section |^ contains results from recent benchmark calculations of the light hadron 
spectrum. In section ^ recent determinations of the strange quark mass are described. 
Finally, section ^ contains some concluding remarks. 

2. Basic concepts 

The lattice formulation of QCD replaces the familiar continuous Minkowski space-time 
by a discretized, Euclidean version with finite volume ■ T . Points in space-time are 
separated by a finite distance a, the lattice spacing. The inverse lattice spacing, a~^, 
acts as an UV cutoff, which regularizes the infinities that are typically encountered 
in Quantum Field Theory. The quark and antiquark fields ili{x), ip^x) are associated 
with the sites of the lattice, whereas the gauge field is represented by the so-called 
link variable [/^(x), which connects neighbouring lattice sites, and is an element of 
the gauge group SU(3). After choosing a suitable, gauge-invariant discretization of 
the QCD action 

S[U, V^, ^] = Sg [U] + Sf [U, ^, V] , (2) 

where Sq is the lattice action of pure Yang-Mills theory and Sp denotes the lattice 
fermion action, one may define the expectation value of an observable Q as 

(O) = I y D[U]D[{p]D[ij] n exp {-Sg - Sp) . (3) 

Here the normalization of the functional integral Z is determined by requiring 
(1) = 1. After performing the integration over fermionic fields the expression for 
the expectation value becomes 

{^) = ^ fli dUf^ix) H dct {D + nif) n exp(-5G), (4) 

where D is the lattice Dirac operator and nif is the mass of quark flavour /. Thus, the 
discretization procedure has given a meaning to the functional integration over gauge 
fields, which reduces to an ordinary multiple-dimensional integral over group elements. 
The various steps leading to eq. (Q) can be treated in a gauge-invariant manner, and 
hence the lattice formulation represents a regularization procedure that preserves the 
gauge invariance of QCD. Furthermore, the definition of physical observables does not 
rely on perturbation theory. Therefore eq. (jij) forms the basis for a non-perturbative, 
stochastic evaluation of expectation values of physical observables using numerical 
simulations. 

Despite an enormous increase in computer power, there remain several major 
difficulties that make realistic simulations of QCD a hard task. Perhaps the 
biggest challenge is the inclusion of dynamical quark effects: the evaluation of the 
fermionic determinant in eq. (Q) in numerical simulations is still very costly, even 
on today's massively parallel computers. In early simulations of lattice QCD the 
determinant was therefore set equal to 1, a choice whose physical interpretation 
corresponds to neglecting quark loops in the evaluation of (ft). Although this 
represents a rather drastic assumption about the influence of quark-induced quantum 
effects, the quenched approximation works surprisingly well, as we shall see later. 
However, it is clearly desirable to develop simulation algorithms that allow for a 
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more efficient evaluation of the fermionic determinant. With present algorithms, the 
cost of simulating "full" QCD is roughly a thousand times higher than that of the 
corresponding quenched simulation. 

An indirect consequence of using the quenched approximation is the so-called scale 
ambiguity. That is, the calibration of the lattice spacing in physical units, [MeV], 
depends on the quantity Q, which is used to set the scale 

a-MMeV] = ^^^, Q = U,mN,mp,... (5) 

This ambiguity arises because different quantities Q are affected by quark loops in 
different ways. Clearly, dimensionful quantities such as particle masses are directly 
affected by the scale ambiguity. 

Another problem one has to address are lattice artefacts (cutoff effects). Let fJ 
denote the mass of a hadron in units of, say, the nucleon mass. Then the expectation 
values on the lattice and in the continuum differ by corrections of order a^: 

{ny^' = (f])™"* + 0(aP), (6) 

where the value of integer p in the correction term depends on the chosen discretization 
of the QCD action. Values of a that can currently be simulated lie in the range 
a « 0.2-0.05 fm. The size of the correction term can in some cases be as large as 
20%, depending on the quantity and the chosen discretization. An extrapolation to 
the continuum limit, a ^ 0, is then required to obtain the desired result. Surely this 
extrapolation is much better controlled if the discretization avoids small values of p. 

Finally there are restrictions on the quark masses ruf that can be simulated. In 
general the following inequalities should be satisfied in any simulation: 

a <C ^ < L, (7) 

where L is the spatial extent of the lattice volume. The quantity ^ denotes the 
correlation length of a typical hadronic state and serves as a measure of the quark 
mass. The inequality on the right places restrictions on the light quark masses that 
can be simulated: if those are too light one may suffer from finite-size effects, since ^ 
becomes large. Typical spatial extensions of L « 1.5-3 fm imply that the physical 
pion mass cannot be reached. The left inequality restricts the masses of heavy quarks. 
Since w 2-4 GeV, it is clear that relativistic &-quarks cannot be simulated. One 
therefore relies on extrapolations in mi to connect to the physical u, d and b quarks. 

3. Hadron spectroscopy 

In many ways the calculation of the mass spectrum of light hadrons is a benchmark of 
lattice QCD. A comparison of the experimentally observed spectrum with the results 
from an "ideal" lattice simulation of QCD, in which dynamical quark effects, lattice 
artefacts and quark mass dependences are all sufficiently controlled, would represent a 
stringent test of QCD as the underlying theory of the strong interaction. Alternatively, 
such a comparison enables us to assess the inherent systematics of current simulations, 
notably the effects of dynamical quarks. In particular one may investigate the quality 
of the quenched approximation. 
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Figure 1. The quenched hadron spectrum from rcf. [[7| compared with experiment 
(dashed hnes). 



3.1. Quenched light hadron spectrum 

Recently the CP-PACS Collaboration presented a precision calculation of the light 
hadron spectrum in quenched QCD Q], which superseded an earlier study by GFll 
The findings of CP-PACS are summarized in the plot shown in Fig. |l|. The main 
features of the spectrum are well reproduced by the quenched lattice data. Bearing in 
mind that only two input quantities have been used, namely the mass of the p-meson 
to set the scale and the mass of a strange meson (either mx or m^) to fix the mass 
of the strange quark, the fact that the masses of so many hadrons are predicted quite 
accurately represents a major achievement in the understanding of QCD. Nevertheless, 
one finds small but significant deviations from the experimentally observed spectrum. 
For instance, the ratio of the nucleon and p masses is calculated as 

mN/mp = 1.143 ±0.033, (8) 

which is 6.7% (2.5cr) below the experimental value of 1.218. Similarly, vector- 
pseudoscalar mass splittings such as tuk* — tok are too small by 10-16% (4-6tT), 
depending on whether mK or is used to fix the strange quark mass. This implies 
that, for the first time, a significant deviation between the quenched QCD spectrum 
and nature is detected. Thus the conclusion is that quenched QCD describes the 
light hadron spectrum at the level of 10%. However, it also shows that the quenched 
approximation works surprisingly well, since the discrepancy is fairly mild. This has 
important consequences for lattice predictions of some phenomenologically interesting 
quantities, for which one has to rely on the quenched approximation for some time. 

Although the CP-PACS results represent a real benchmark in terms of statistics, 
parameter values and lattice volumes, further corroboration of these findings 
is required. Recent calculations employing different discretizations have largely 
confirmed the findings of Q: the MILC Collaboration Q has used staggered 
fermions pO{ and finds a value for mN /mp in the continuum limit of 

niN/mp = 1.254 ± 0.018 (stat) ±0.028 (syst). (9) 
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This is in broad agreement with experiment, but the difference to the CP-PACS 
result amounts to only 2a. A recent calculation by UKQCD using 0(a) improved 
Wilson fermions |11| finds m^/nip = 1.2Q'^^^. Whereas the error is too large to 
detect a significant deviation, UKQCD 's results for the spectrum also indicate that 
the quenched light hadron spectrum agrees with experiment at the level of 10%. 



3.2. Beyond the quenched approximation 

An obvious question is whether sea quark effects can account for the observed 
deviation of the quenched light hadron spectrum from experiment. Recently several 
collaborations have studied the light hadron spectrum using N{ = 2 flavours of 
dynamical quarks, which are identified with the physical up and down quarks ]T^-p^. 
However, the masses of the dynamical quarks that can be simulated are quite large. 
A measure for the difference between the light quark masses found in nature and 
those used in dynamical simulations is provided by the pseudoscalar-to-vector mass 
ratio, TOps/mv. Current simulations are typically performed for mps/my = 0.6- 
0.8, whereas the physical value is mpg/my — mj^/mp — 0.169. Therefore one relies 
on extrapolations in the sea quark mass to make contact with the physical case. In 
addition, it has not been possible to simulate a third dynamical quark with mt « mg 
efficiently using known algorithms. Therefore, in most dynamical calculations an 
unphysical number of sea quarks has been used. 

Despite these shortcomings one can perform systematic studies of the influence 
of dynamical quarks on the light hadron spectrum. CP-PACS ||l^ have studied the 
continuum limit of hadron masses computed for iVf — 2 and compared it with the 
results of the quenched light hadron spectrum discussed before. They find that the 
discrepancy with experiment is reduced from 6.5% to 1.4% for computed with 
N[ — 2. For tok» the gap decreases from 4.4% to 1.0%. This demonstrates clearly 
that the experimentally observed spectrum is reproduced more closely when sea quarks 
are "switched on" . Furthermore, results for the vector-pseudoscalar mass splitting 
reported by UKQCD ||l^ also show that lattice data for this quantity approach the 
experimental value as the sea quark mass is decreased. 

Of course, the small remaining differences between QCD with iVf = 2 and 
experiment have to be explained. It is reasonable to assume that the unphysical value 
of N{ will have some influence. Further studies are also required to decide whether 
there is yet sufficient control over the extrapolations in the quark masses and those to 
the continuum limit. 



4. The mass of the strange quark 

Quark masses are important input parameters in many theoretical applications 
but, despite a great deal of activity, relatively little is known about their absolute 
values ijl^ . Here I shall describe calculations of current quark masses defined through 
the PCAC relation. For charged kaons it can be written as 

/Km| = (mu + ms)(0|u75s|K), (10) 

where the bars above the quark masses indicate that the "running" masses are 
considered, which depend on the energy scale. In order to determine the sum of 
up and strange quark masses using the experimental result for /kw^ one simply has 
to compute the matrix element (0|u75s|K) in a lattice simulation. However, in order 



Strangeness in Lattice QCD 



6 



Lattice 


non-perturb. 


scheme: miat(a) 





Intermediate 
scheme: tox(a*) 
Mo < M < M 



perturbative 



MS scheme: 



"MS 



(m) 



Figure 2. Schematic relation between quark masses in lattice regularization and 
the MS scheme through an intermediate renormalization scheme X. 



to represent meaningful theoretical input, quark masses have to be renormalized, and 
hence their values depend on the adopted renormalization procedure. By convention 
quark masses are quoted in the MS scheme of dimensional regularization at a reference 
scale, typically /i = 2 GeV. This implies that the renormalized matrix element in the 
MS scheme has to be determined; it is related to its lattice counterpart by 

{0\u^5s{p)\K)^ = zf^{go,afl){0\u^,s\K),,f (11) 

Here go is the bare coupling, and fl is the subtraction point in the MS scheme. The 
renormalization factor is known to one-loop order in lattice perturbation theory. 
However, it is well known that lattice perturbation theory converges slowly, and in 
order to remove all doubts about the reliability of the matching procedure it is evident 
that a non-perturbative determination of the renormalization factor is required. 

However, a non-perturbative normalization condition relating the matrix elements 
in lattice regularization and the MS scheme cannot be formulated, since the latter is 
only defined to any given order of perturbation theory. This technical difficulty can 
be overcome by introducing an intermediate scheme X and considering a two-step 
matching procedure, as shown schematically in Fig. |^. The first part is the matching 
of the bare current quark mass miat(a) to the running mass in the intermediate scheme 
rHxip)- This amounts to computing Zp ((70,1/^0) between the lattice scheme and 
scheme X at a fixed scale /io- The second part is the determination of the scale 
dependence of the running mass tox(m) from /xq up to very high energies, where the 
perturbative relation between tox and is expected to be reliable. Through this 
two-step process the use of lattice perturbation theory is completely avoided. 

Two examples for such intermediate schemes have been proposed. The so-called 
"regularization-independent" (RI) scheme is described in Another intermediate 
scheme is defined using the Schrodinger Functional (SF) of QCD This scheme 

allows one to compute the scale dependence non-perturbatively over several orders of 
magnitude, using a recursive finite-size scaling technique. Once the scale dependence 
of the running mass in the SF scheme, msF, is known up to /i « 100 GeV one can 
continue the scale evolution to infinite energy using the perturbative renormalization 
group functions and thereby extract the renormalization group invariant (RGI) quark 
mass M. At this point the matching to the MS scheme is trivial, since M is scheme- 
independent, i.e. it coincides in the SF and MS schemes. 

Figure ^ shows the non-perturbatively determined scale dependence of rfisF/M 
computed in quenched QCD and compares it to the perturbative scale evolution. 
The left-most data point in Fig. ^ corresponds to a scale /xq — 275 MeV. It is 
here that the matching between lattice regularization and the MS scheme via the 
SF is completed, through a non-perturbative calculation of Zp^{gQ,a^o) at fixed 
fio = 275 MeV. All dependence on the intermediate SF scheme and the scale fio 
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Figure 3. Non-perturbative scale evolution of mg-p/M computed in lattice 
simulations of the SF (solid circles). The lines correspond to the scale evolution 
computing using various orders of perturbation theory. 



drops out in the total renormalization factor 
wsf(mo) 



M 



(12) 



where msF{fJ-o)/M — 0.864±0.011, as can be read off the left-most pomt in Fig. || [ 
The total renormalization factor is thus known with a precision of 1.5%. 



The strange quark mass in quenched QCD 

I shall now discuss the application of non-perturbative quark mass renormalization 
to compute the mass of the strange quark in the quenched approximation. Before 
discussing the details, it is useful to recall that mass ratios of light quarks are predicted 
quite accurately by Chiral Perturbation Theory (ChPT), for instance |2^] 

^ = 24.4 ±1.4, M = i(Mu + Afd). (13) 

In order to determine Mg, it is then sufficient to compute the sum {Mg + M) on the 
lattice and combine it with eq. (p^). 

Here I shall concentrate on the calculation performed by the ALPHA/UKQCD 
Collaboration |2^. By combining the PCAC relation, eq. (p^), with the total 
renormalization factor, one obtains the sum of RGI quark masses (Mg + AI) in units 
of the kaon decay constant| 

/k V^SF ZpJ (0K75s|K)iat 

The dependence on the lattice spacing can be eliminated by performing a continuum 
extrapolation as shown in Fig. ^ Using the experimental value /k = 160 ± 2 MeV Jl7[ | 
one obtains the result for the sum of RGI quark masses in the continuum limit: 

Ms M = 140 ± 5 MeV. (15) 

§ Here, £ denotes the Dirac spinor for the isospin-symmetric light quark. 
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Table 1. Estimates for the quark masses rris and m in the MS scheme in 
quenched QCD. The choice of intermediate renormahzation scheme X is also 
shown, where applicable. 



Collaboration 

CP-PACS [0 
Becirevic et al 
Gockeler et al 
Wingate et al 
ALPHA/UK 
Blum et al |29| 
JLQCD [0 
Becirevic et al jplj 
Gimenez et al 132[ 




m{2GeV) 


ma(2GeV) 


X 


a — > 


4.4(2) 


110(4) 




yes 


4.8(5) 


111(9) 


RI 


no 


4.4(2) 


105(4) 


SF 


yes 




130(11)(18) 


RI 


no 




97(4) 


SF 


yes 




96(26) 




yes 


4.23(29) 


106(7) 


RI 


yes 


4.5(4) 


111(12) 


RI 


no 


5.7(1)(8) 


130(2)(18) 


RI 


no 







This result can now be converted into m^^(/i) at /I = 2 GeV. First, eg. ([Tsl) is 
combined with the prediction from ChPT, eq. (|l3|). Then, the relation between the 
running mass in the MS scheme and the RGI quark mass is computed in 4-loop 
perturbation theory as 

m^(/i)/Af = 0.7208 at /i = 2GeV. (16) 



This procedure yields the final result in the quenched approximation 




mf^(2GeV) = 97±4MeV. (17) 

The quoted uncertainty of ±4 MeV contains all errors, except those due to quenching. 
The high precision of this result is a direct consequence of recent progress in lattice 
calculations, in particular the implementation of non-perturbative renormahzation and 
control over lattice artefacts. As mentioned in section |^, the conversion into physical 
units is ambiguous in the quenched approximation. For m^^(2GeV) the resulting 
uncertainty was estimated in ref. ||2^ to amount to ~ 10%. 

Table |l| contains a compilation of recent results for TOs and fh obtained in the 
quenched approximation. Direct comparisons of these results should be made with 
care, since systematic errors have not been estimated in a uniform manner. Also, 
the conversion into physical units has been performed using different quantities. This 
manifests itself in the typical spread of the central values, which is of the order of 10%, 
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consistent with the above estimate. It is still quite remarkable that lattice estimates for 
light quark masses have stabilized, which represents a big improvement with respect 
to the situation before 1998. 

4-2. Sea quark effects in rus and fh 

So far the most comprehensive study of quark masses using simulations with dynamical 
quarks has been presented recently by CP-PACS Ea]. Earlier results can be found 



in refs. \12, 3^. When discussing the results, it is important to keep in mind 
that non-perturbative renormalization has not yet been implemented for N{ = 2 
flavours. Instead the renormalization of the lattice matrix element is usually performed 
using so-called "mean- field improved" lattice perturbation theory |^^. In addition 
to computing the strange quark mass, CP-PACS have also estimated the average 
light quark mass through extrapolations to the chiral limit. Furthermore they have 
compared different definitions of the quark mass, based either on the axial vector Ward 
identity (AWI, i.e. the PCAC relation) or the vector Ward identity (VWI). 

The continuum extrapolations for the strange quark mass are shown in Fig. ^. 
Results from the AWI and VWI definitions of the quark mass have been extrapolated 
enforcing a common continuum limit. By comparing the results for A'f = 2 to those 
obtained in the quenched approximation, one finds that the inclusion of dynamical 
quark effects decreases the estimate for the strange quark mass by roughly 20%. The 
results for nis by CP-PACS can be summarized as follows: 

The light quark mass fh shows a similar decrease when sea quark effects are included. 
However, for both Ni = 2 and iVf = (quenched) one finds that rris/fh « 25 in the 
continuum limit, which is consistent with the prediction from ChPT, eq. (p^). 

Figure |^ also illustrates that lattice artefacts are quite substantial, so that an 
estimation of quark masses at non-zero lattice spacing would be misleading. This 
underlines once more the importance of the continuum extrapolation. 
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5. Concluding remarks 

Lattice simulations of quenched QCD have reached a level of precision of a few per 
cent, thanks to increases in computer power and conceptual advances. The picture that 
has emerged is that the quenched approximation works surprisingly well, as signified 
by the success in predicting the spectrum of low-lying hadrons. Non-perturbative 
renormalization and good control over lattice artefacts are the crucial ingredients that 
lead to precise estimates of the mass of the strange quark. Also, in order to control 
the region of very light quarks, which cannot be simulated on present machines, it 
is useful to combine lattice simulations with Chiral Perturbation Theory. In fact, a 
combination of these techniques may be instrumental in resolving the long-standing 
question of whether the up-quark is massless . A pilot study along these lines has 
already appeared in the literature |3^]. 

A major share of computational resources is devoted to studying the effects 
of dynamical quarks. First results indicate that sea quarks "are doing the right 
thing" , since the discrepancies in the light hadron spectrum between simulation and 
experiment are decreased. However, more effort is required to understand sea quark 
effects on a quantitative level. This not only requires larger computers but also the 
implementation of additional technology such as non-perturbative renormalization for 
the dynamical case. Finally, there are algorithmic challenges: the efficient simulation 
of sea quarks for which mps/my <C 0.5, and also of odd Ni is crucial for the 
development of lattice QCD into a fully fledged phenomenological tool. 
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